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(54) Abstract Title 

Signal arrival direction estimation . 

(57) In a DS-CDMA communication system for performing a transmission by a spreading process using a 
code having a period longer than a . symbol length, received signals from an array antenna 1-0 to 1-{M-1) are 
subjected to a despreading process by despreading units 3-0 to 3-(M-T>. A cross-correlation function with time 
lags of the received signals is obtained by Inverse modulating despread output signals by a known symbol In 
an inverse modulator 8. An arrival direction of a received signal is estimated by an arrival direction estimation 
unit 6 based on the cross-correlation function. A beam former 7 combines the despread output signals 
depending on the arrival direction and outputs an output signal to a channel receiver The arrangement 
provides an accurate estimate of the arrival direction of a desired signal regardless of reception level and high 
interference. 
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At least one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy. 
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ARRIVAL DIRECTION ESTIMATION METHOD USING AN ARRAY 
ANTENNA AND PS -CDMA RECEIVER UNIT USING THE METHOD 

5 " " : ~ •• • • " ~ ' 

The present invention relates to an arrival direction 
estimation method using an array antenna for estimating an 
10 arrival direction of a desired signal. The present 

invention further relates to a DS-CDMA receiver unit in a . 
Direct Sequence Code Division Multiple Access (DS-CDMA) 
system. 

In a base station of a DS-CDMA communication system, 
15 signal, reception is performed by utilizing . an array antenna. 
Further, an arrival direction of a desired signal is 
estimated from signals received from the array antenna. By 
. . transmitting and receiving in a beam directed in^ the. - 

estimated arrival direction, it is possible to reduce both 
20 interference and required transmission power due to an < 

increase in antenna gain. Therefore, it is desirable -that 
the arrival direction of the desired signal is estimated 
. accurately. . V; 

F^J« 7 of the a:ro ipsrying dradngs is a diaLydn Viewing a receiver unit cf a teee 
25 station using an array antenna. In particular/ Fig. 7 shows 
array antenna elements 71-0 through 71- (M- 1 ) forming a-' . 
/ receiving antenna, receptibn processors 72-0 to 72- (K-l) , 
despreading units 73-0 to 73- (M^l) > beam forming unit 74, a 



channel receiver 75, an arrival direction estimation unit 76 
and a beam former 77. 

Received signals r ,0> to r^ 1 3 from the array antenna 
elements 71-0 to 71-(M-1) are input to the reception 
5 processors 72-0 through 72- (K-1) corresponding to the 

different channels. In addition, the reception processor 
72-0 carries out. a despreading process with respect to the 
received signals by . the despreading units 73-0 to 73-(M-l), 
which are input with a spreading code Co. The despreading 

10 units 73-0 to 73-(M-l) provide despread output signals x ,0} 
to x 1 * 1 " 1 * to the arrival direction, estimation unit 7 6 and. the 
beam former 77 of the beam forming unit 74. The arrival, 
direction estimation unit 7 6 estimates an arrival direction 
of a desired signal based on a cross-correlation function 

15 among the received signals received by adjacent antenna 
' elements . 

The beam former 77. multiplies a. weighting coefficient 
obtained by estimating the arrival direction by the arrival 
direction estimation unit 7 6 to the despreading -output 
20 signals x <0> . through. x (M "^ \ i Further;, the beam former 77 then 
combines . the products to provide an output signal . to the 
channel .receiver, 7.5. Therefore, received data of the ; 3ctti 
, channel is re-generated .'by the reception processors 72r0 to. 
72- (K-l) corresponding to channels' 0 through K-l ; 
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Fig. 8 of the acrxiTperying dradrj^ is a diagran shewing a pcejicxjs3y<xxisidaoed b^an 
forming unit 74 including the arrival direction estiination unit 76. 
In particular, Fig.. 8 shows the beam forming unit 74 
including an arrival direction estimation unit 76 and a beam 

.5 former 77 having. multipliers 78-0 to 78-3 and an adder 79. 
Received signals x <0) (n) to x f3> (n) from the antenna 
elements 71-0 to ?l-3 (refer to Fig. 1) are input to the 
arrival direction estimation unit 76 and the beam former 77. 
The arrival direction estimation unit 7 6 outputs to the beam 

10 former 77 weighting coefficients aO to a3 obtained by the 
estimation of the arrival direction of the desired signal. 
Th^ weighting coefficients aO to a3 are respectively input 
to the multipliers 78-0 to 78-3 of the beam former 77 to be 
multiplied with signals x f0> (n) to x (3} (n) and then coherently 

15 summed by the . adder 79. An output signal y(n) of the adder 
79 is then input to the channel receiver 75 (refer to Fig. 
7) . 

Fig. 9 of the auua ip a r y ing dcaringg is a diagram stewing a pt^icxriy-cx^nsddecBd arrival 
direction estimation unit . In particular/ Fig. 9 shows an 
20 arrival dlre^ unit 76 including correlation 

calculators 81-6. through 81-2, an adder 82/ and an array 
weighting, coefficient calculator 83. A correlation . 
calculator 81 is shown in the iowei? part of Fig. 9, wherein 
• each of the correlation calculators 81-0 to 81-2 include a 
25 ..multiplier 84 and. an averaging f i 1 1 e r 85. Further/ received 
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signals x. (m> (n) and x (re+11 (n) from mth and (m+l)th adjacent 
antenna elements are input to the multiplier 84. In this 
case, one of the received signals is a complex conjugate 
'(indicated by.*), and a complex conjugate product between 
5 'the: received signals from, the adjacent antenna elements 

(cross-correlation function of the zero lag) is obtained, A 
temporal average is also performed by the averaging filter 
85 so tbat a cross-correlation function R (ml is output. 

A received signal r fm} (t) from a mth antenna element 

10 at a time t satisfying (n-I)t < t < nT is described by the 
following formula (1) , where 0 ± {my for 1=0 to N-l in 
formula (1) is described by formula (2) f N denotes a number 
of Users/ Ai . denotes a received amplitude of an ith user 
signal, c A (t) denotes a spreading code of the ith user, Ti 

15 denotes a relative delay of the ith user, h x denotes an ith 
user transmission symbol, d denotes a distance between the 
antenna elements,. 6* denotes an arrival angle of the ith 
.. user sigrnal, . N (tt) <t) denotes a noise signal, T denotes a 
symbol length, and a multipath is not taken into 

20 consideration for the sake of convenience - 

r M (t)« Z AiCitt-TilbCt-TiJexpC-DC^M+N^ (t) 

;> : ;V.v ;^-'^ c ^^27i^sin8t • . • . (2) 
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A despreading output signal x a) (m> (n) which is obtained 
by subjecting the received signal r <B) (t) from the mth 
antenna element to a despreading process by the spreading 
code p, u (t) of the kth user is described by the following 
5 formula (3) , where w± } (n)in the formula (3) is described by 
the formula (4) . 



10 



nT 



Xk 1 "' (n-)« 



C (t)c K * (t-T k )dt 

J (n-l)T 



15 



20 



(nT 



|(n-l)T 



■n-i 

XAi 



i=0 
i^k 



Ci (t-Ti)c k * {t-x k )bi (t-Ti)exp(-j0 tn ") J'dt 
+N (B, .(n) ... 



25 



30 



N-l nT 
=A k bK(n)exp{-j0 {,,,) >4- ZM f 

i=0 J <n-l)T 

i*k 

Ci (t-tj) c k * <t-T k )bi (t-Ti)dt] exp (-jo'") 
+N ,B » (n) 



35 



SUbk (n) exp (- Is?**}* 
"N-l ',. 

Z [Ai w~(n)exp M 0i ,m> ) ] +N <B,) .(n> 
i=0 

i*k ../'. 



— : (3) 
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Wij(n)= \ dtt-XiJcK * (t-T K )bi (t-xjdt 

J (n-l)T 

5 >~ (4) 

In addition > the correlation calculator 81 of Fig.. 9 
calculates a correlation R x tm> between the received symbol 
10 x ,ml (n) of the mth antenna element and the received symbol 
x lB4U (n) of the (m+l)th antenna element bas.ed on the 
following formula (5), where brackets **< >" indicate 
averaging and the terms inside the brackets "< >" are shown 
. in .the formula (6) . 



15 



Rr^'^x^' (n)x k (B,+1, *(n)> — (5) 



Xk*"" Cn)x fc ^ l, *(n) = lA,,| 2 |bK eitpIj-^-Bk^Jj ' 

N-l 

20 / .+ E |A i | 2 |w ik (n)| 2 exp[j{0 lc tm) . 

' ' i«o • ; •• "• ■■" 

i*k 

-0k'"*") )+N ,,B3 *{n)N ,B, *n 

(6) 

25 . ' - : 

The' correlation R k (a " is then averaged for each of the 

.antenna elements and a cross-correlation function. R k among 

the adj acent antenna elements is obtained based on the 

.following formula (7) , where M denotes the number of .. antenna 

30 . elements. 



• .' t : .. ■ .' .' M-2 . . -. 

Rk- [i/ (M-l) ]X Ric 0 " . — (?) 

35 ' 
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Fig. 9 shows a structure for adding the correlations 
R i0) to- R (2) from the correlation calculators 81-0 to 81-2 by 
the adder 82. However/ no calculation function related to 
. the term [1/<M-1)3 in the formula (7) is shown. -Since (M-l) 
5 is a fixed value, it is clear that an averaging process can ; 
be carried put based on the addition. 

In formula (6) as described above, the second term on 
the right hand side of the equal sign is an uncorrelated 
signal and the term {w^fn) I is a small value* Hence, if 
10 there is no significant difference ajriong the received signal 
amplitudes Ai of each of the users, it is possible to 
extract by averaging a phase error of a desired signal (kth 
user) of the first term on the right hand side of the equal 
. sign of the formula (6).. Accordingly, based on a ratio of. , 
15 an imaginary part Im(R k ) and a real part Re (Rr) using the 
cross-correlation function R k of the formula (7), it is 
possible to estimate, from the following formula (8), an 
y arrival angle 8 k of the kth user signal* 

20 e^sin^t (l/27cd)tan" 1 {Im(R k )/Re{Ri t ) 1] — (8) 

According .to the previously-considered arrival direction . 
estimation method described, above, the arrival direction is 
estimated from the strongest signal received after the 

25 -despreading process, regardless ; of whether or not the . 

desired signal exists + For this, reason:/ - it is difficult to 
estimate the correct 'arrival direction of the desired 'signal 



in cases where the reception amplitude of the desired signal 
is small and in cases where large interference signals 
exist. 

.For example, when a multi-rate transmission including a 
5 high-speed rate and a low^ speed rate, in the DS-CDMA 

communication system is made, the received signal of the 
high-speed rate user becomes an interference signal for the 
received signal of the low-speed rate user since the. 
transmission power of the high-speed rate signal is large. 
10 Consequently, a problem exists in that it is difficult to 
estimate the arrival direction of the desired signal of the 
low-speed rate user. 

15 ■ 

Accordingly, it is desirable? ■ : — ~ — : — — — 

to. enable estimation of the correct arrival direction of the 
desired signal in a DS-CDMA communication system using a 
code having : a period ].o.nger than the symbol length of the 

20 spreading .cpcte-/ : 'ey^nf.wH^n:-' a/ ■ i^t er f er'enae _ : : s ighai- exists. 

An embodirrient : of p first aspect — , . • ; . -. — — - — — ■ 

• - — ■ . • ■ - - h -~ of the present 

: invention, uses ; an array antenna in a DS-CDMA communication 
system for performing a transmission by performing a 
25 spreading prod ess including a code having a period longer 
than a symbol length. The method also includes to subject 



received signals from the array antenna to a despreading 

process • In this errtrriimsnt a cicssr<3QfcXTBlartiQn fixcticn of lacp/ p.dLtLLobLy excluding a aeto 

time lag, for signals corresponding to adjacent antenna 

elements is obtained by inverse-modulating despread output 
5 signals by a known symbol- Further, an arrival direction of 
a desired signal based on the cross-correlation function is 
estimated. 

In ax>t±er ertrrliment it is alpn prRPjihlp to use ava^rrrp valine cf the 

cross-correlation function for a plurality of different lags 

10 as the cr o ss - corre lat ion function of the lags among the 

inverse modulated signals corresponding to the adjacent 

antenna elements, . wherein the zero jag is excluded from this 

function* It is also possible to obtain the cross- 

: correlation function using both positive .and negative lags 

15 excluding the zero lag. among the inverse modulated signals 

corresponding to the adjacent antenna elements. A pilot 

symbol may be used as the known symbol to iiiverse-modulate 

the received signals ? It is also possible to provide a 

decision feedback symbol as the known syiabol to inverse- 

20 modulate the received signals. 

An effbtVii merit of , a second -aspect of the present invention can, provide' 

; • : a^ re-03MA receiver unit for use in. a DS-CEMA' ccmnunication 
• ;. iYstem,v which: i carries but a transmission by performing a 
' ;: -\\ spreading, process using a. code having a period longer than , a 
25 syrtbl Isngtiv In this grtacdinent the Z&<n&xazewer unit irciudes de^esdiig units foe 



9 



'-''-subjecting received signals from array antenna elements to a 
despreading process. Furthar included in this enfoodiment is a beam forming 

unit 4 .for inputting despread output signals from the 
■despreading units and a channel receiver for inputting •* 
5 output signals of the beam forming unit. The beam forming 
unit includes an inverse modulator for inverse-modulating 
the despread output signals by a known symbol/ an arrival .' • 
direction estimation unit for receiving signals from the 
inverse, modulator and a beam former. The beam former 
10 provides an output by multiplying a weighting coefficient 
from. the arrival direction estimation unit to the despread 
output signals and then adding such products. 

The inverse modulator may be configured to. use either a 
pilot symbol or a decision feedback signal, as the known 
15 symbol. The arrival direction estimation unit may include a 
correlation calculator for obtaining a cross-correlation 
. among output signals of the inverse modulator corresponding 
. - to the "adjacent antenna elements and an array weighting 
coefficient calculator for obtaining and averaging the 
20 cross^correlation functions with respect to a plurality of 
. lags excluding a zero lag. The array weighting coefficient 
calculator then outputs a weighting coefficient obtained as 
a result of estimating the arrival direction based thereon. 
In this case, the cross-correlation function with respect to 
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the plurality of lags may be calculated by providing first 

and second shift registers. 

Reference will now be made f by way of example only, to 

the accompanying drawings, in which: 

Fig.. 1 is a diagram showing ~ ■ — a DS-CDMA 

receiver unit embodying the present invention; 

Fig. 2 is a diagram showing a spreading modulator of a 

transmission unit for use with an enrbodiment of the present invention; 

10 Fig- 3 is a diagram showing a despreading unit 

for use in an enribcdiment of the present invention; 

"Fig;.". 4 is a diagram. showing a beam forming unit 
for use in an errbodirnent of the present invention; 

Fig. 5 is a diagram showing an arrival direction 
15 estimation unit for use in an embodiment of the present invention; 

Fig » 6 is a diagram showing a correlation calculator 
for use in an errbodLnent of the present invention; 

Fig. 7, discus9ed hereirbefbrei, is a diagcan shewing a ledeiver init of a beee 
station using an. array antenna; 
20 Fig. 8, < fti?gT3=g=d hereirbefege, is a diagran shewing a pgevin Kly-<TTHidfingd bean fenning 

• . unit haying estimation unit; "and 

Fig.: 9, <3dsb yyyri hereiitefEge, is a diagram showing a preyj a gly-^u tiidggdl arrival 

• /' % direction e'stimation -unit. • : ; 
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Fig. 1 is a diagram showing an embodiment of a DS-CEMA 
receiver unit according to the present invention. As can be 
5 seen from Fig- 1, the receiver unit includes M array antenna . 
elements 1-0 to 1- {M-l ) forming a reception antenna, K 
reception processors 2-0 to 2-(K-l) provided in 
correspondence with K channels, M despreading units 3-0 to 
3t(>4^1), a beam forming unit 4 and a channel receiver 5. 
10 The beam forming unit 4 includes an arrival direction 

estimation unit 6, a beam former 7, and an inverse modulator \- 
-•; 8. The reception processors 2-0 to 2- (K-l) each have the 
■same structure and receive signals from the array . antenna 
elements 1--0 to. 1.- (M-l) • 
;.. ' is Fig. 1 further shows the receiver unit of the DS-CDMA 

communication system using a code C 0 having a period. longer 
than a symbol length. Received signals r (D1 to r CM ~ 1} from the 
M array antenna elements 1-0 to 1-fM-l) are respectively, 
input, to the despreading units 3-0 to 3 (M-l) in the 
V; 20 reception processor 2-0 . For example, the despreading is _ : 
carried out by the use of a spreading code C 0 . 

; In this case, a' transmission unit, of -a mobile Station 
:/V : or the like would include a spreading modulator for the . 

- code. Fig. 2 shows an example of such a spreading . 
. 25 ^bduiator/ which includes multipliers {exclusive-OK . . 
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circuits) 11, 12 and 14 to 17, a scramble modulator 13, 
polarity converters. 18 to 21, and adders 22 and 23. 

With respect to a data symbol, a short code SCi of the 
quadrature code is multiplied by the multiplier 11 to 
5 provide code spreading. With respect to a pilot symbol, a 
short code SC q of the~ quadrature code is multiplied by the 
multiplier 12 to provide code spreading. As a result, a 
distribution is made to quadrature channels (I,Q). 

In addition, codes LGi and LC q of the quadrature code 
10 which have a period longer than the symbol length are also 
distributed to the quadrature channels (I, Q) and input to 
the multipliers 14 through 17. For example with respect to 
the data symbol which is subjected to code spreading by . the 
short code SCi, the multiplier 14 multiplies .the code LCi 
15 and the multiplier 17 multiples the long code LC q to carry 
. out a scramble modulation. With respect to the pilot symbol 
which is subjected to code spreading by the short code SC q/ 
the multiplier 15 multiplies the code LCq and the multiplier 
16 multiplies the long : LC 4 to also carry put scramble 
20 modulation . In this case, the pilot symbol may be inserted 
and transmitted for every predetermined length of data 
symbols, as a known symbol having a number of symbols less 
than a predetermined length. Alternatively, it is possible 
ytp transmit the pilot symbol.. together with the data symbol 
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■' having a phase that is orthogonal with respect to the data 
symbol 

Output signals of each of the multipliers 14 to 17 are 
converted to the polarity converters 18 to 21. For example, 
5 a logic n O n .is converted into "+1" or a logic n l n is 

converted into "-1% which are then input to the adders 22 
and 23.. The adder 22 adds an output signal of the polarity 
converter 18 with a polarity, and adds an output signal 
of the polarity converter 19 with a polarity. In other 
10 words, the adder 22 carries out a subtraction process . In 
addition/ the adder 23 sums the output signals of the 
. • polarity converters 20 and 21. Accordingly, output signals 
of the adders 22 and 23 take a ternary value of +1, 0 or -1. 
The ternary value signals from the adders 22 and 23 are 
15 transmitted after carrying out a quadrature modulation or 
the like in a radio transmitter, which is not shown in Fig. 
2. The radio transmitter transmits signals subjected to 
spreading modulation by the code. 

In the embodiment of the DS-CEMA .receiver unit shown in 
20 Fig. 1, the despreading units 3-0 to 3(M-1) carry out 

;p the process . of the spreading- 

modulatory shown ; in ; Fig.. 2 . As shown .in Fig; 3, for example, 
... each despreading unit includes a descr ambler demodulator 31 ' 
'.. and despreading units 32 tp 35 . The . descr amble demodulator 



14 



31 includes multipliers 36-1 to 36-4, adders 37-1 and 37-2. 
As can be seen from Fig. 3, all of the inputs and 
, ' outputs include 1 and Q channel signals. However, in Fig. 

.1, such detail was omitted. Instead, in Fig. 1, each of the 
5 I and Q channel, signals are represented by a single signal. 
For example, receive signal r{0), code C 0/ "data signal D(0) 
and pilot signal x(0) of Fig. 1 correspond to I and Q 
■ channel signals, codes Ld and LC qr data signals Di and D Q , 
and pilot signals xi and x Q , respectively in Fig. 3. 
10 The received I and Q channel signals are first 

subjected to the quadrature modulation in the radio 
receiver, the illustration of which is omitted. The same I 
and Q channel signals are then input to the descramble 
demodulator 31 to be descrambled using the codes LCi and 
15 LC q . The descrambled signals are then input to the 

despreading units 32 to 35. Two of the descrambled signals 
are subjected to despreading using the short code SCi in the 
despreading; units 32 and 33, which respectively output data 
symbols Di: and D Q of the I and Q channels. While the other 
V-2Q*-; -i- fcwb descrairibled signals are subjected to despreading using 

the short -code SC q in the despreading units. 34 and 35, which 
es^ecti^elyvput^ut •pilot ; teyi^oisi -xt -arid, xq of the I and Q - 
* channels.-.-'. -■ \" . • • • 
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As described above, the despreading units 3-0 to 3- (M- 
1) of Fig.. 1 are configured to output the data symbols and 
the. pilot symbols, by carrying out the despreading process* . 
In the present invention, the cross-correlation function 
5 which is used to estimate the arrival direction is obtained 
using the known symbol. This known symbol may be the pilot 
symbol, as described above. Alternatively, it is also 
■/■possible to feed back a data symbol to the inverse modulator 
8 and to use this data symbol as the known symbol. In 
TO feeding back this data symbol, it would be judged by a 

judging: unit included in the channel receiver. Fig. 1 shows 
•the case, where the pilot symbol is used, as the. known symbol . 

.. . When signal x k (m) (n) subjected to the de spreading is : 
further subjected to inverse modulation using a known 
15 transmission symbol b*(n) of the kth user, a symbol is 
denoted by z k tn] (n) , which is described by the following 
formula (9), The formula (9) is developed into the formulas 
.- (10) and (11} / where W ±j (n) in formula (11) is defined by 
the following formula ( 12 ) . 



20 



\ 25 



z,c cm) (h) =x k (ro) (n)bk* <n> — - (9) 

' V\; "' jit • N-i 

z^ m, ;(n)=A k exp(-jcr fta) )+ . C [TAi*(M0. 

J (n-l)T i=0 

:. : ; : . i*k • . 

: ; c H * ( t-xx) exp {- j 0 tm> ) ] dt+N^ (n) 
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N-l JlT 

=A Jt exp(-j0 1 - ) )+ ZAi[ ( c<t-Xi) 
i»0 ' (n-l) T 



-—(10) 



N-l 

Zk^' (n)aA k exp (- ja*"' )+ 2 [A lWl3 (n) expf-jV 8 ') 

i=0 
i*k 

+N'"» (n) — (ii) 

■ hT 

Ci (t-T t )c k * (t-T k ).dt 

— (12) 



wij(n)= \ Ci (t-r t ) Ck * (t-* k )dt 

J (n-l) T 



in formulas (9) through (12), A* denotes a reception 
amplitude of the kth user, c k (t) denotes a spreading code of 
the kth user, T k and x x respectively denote relative delays 
of the kth and ith- users, N <B " (n) denotes a noise signal, and" 
T: denotes a symbol length. 

.As. may be seen from the following formula (13) , a 
correlation R k ,m ? between a reception symbol z ( "" (n) '• . 
corresponding to the mth antenna element and L symbols 
before and" after the reception symbol corresponding to the. - 
(m+1) th antenna ' element, is obtained by multiplying (1/2L) a 
summation of complex conjugate products amounting to 2L 
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symbols in the positive and negative directions from -L to 
+L, excluding the 2ero lag, and averaging. The terms 
following the summation Z are shown in formula (14) . 

5 R k to, =<(l/2L) S z M (n) z' 6 - 11 * (n+l)> 

1=-L 

1*0 ' (13) 

iO. Z z <»> (n) a 1 "*"* (n+1) 

- 1=L 
1*0 

L N-l 

15 =[A k | i exp [j(0x (o,) -0K (n * ,, ))+ 2 Z|Ai| 2 W ik (n) 

1=-L i=0 

1*0 i*k . 

20 W tk * (n+l)exp[j (e> 1 < " ) -0 i <,,,+^, ) ] 

L 

+ 2 N <B " (n)N |B> * (n+1) (14) 
1=-L 
25 1*0 

The formulas (13) and (14). show a case where both the 
positive and negative lags of the lags (time delays) other 
: than the zero lag are used, where the correlation is 
30 obtained for the L symbols with respect to each qf the: 
. positive and negative lags. When obtaining the cross- ' 
correlation for a plurality of different lags (time delays), 
it. ".is possible to obtain the correlation for the L symbols 
of a plurality of lags in only the positive direction or . 
35 negative direction. In addition, it- is also possible to- 



obtain the correlation for different numbers of symbols of 
different lags in the positive and negative directions • 

As described above, the correlation is obtained for L 
- symbols of the different lags in each of the positive and 
5 negative directions , excluding the zero lag (time delay), of 
the reception symbol z (al (n) corresponding to the mth antenna 
element and the reception symbol 2 <to4l) (n) corresponding to 
the (m+l)th antenna element- The arrival direction is 
estimated from this correlation. Hence, the_ arrival 
10 direction can be estimated by suppressing undesirable 

effects/ even in cases where a phase deviation occurs due to 
" fading. 

By obtaining the average of R <m1 by formula (13), a 
...cross-correlation function among the adjacent antenna 
15 elements can be obtained from the following formula (15). 
'■' In addition, an arrival angle 8 k of the kth user can be 

obtained from the following formula (16) which is similar to 
formula (8) as described above, using the cross-correlation 
. function Rk- ': /■'.;*. . ' 

:r^{i/(m-i)} Z :R k w . — - (is) 

■ : p^o .. 

25 ' ■ flj^ain" 1 [ ll/27id) tan' 1 { im (R^) /Re (Rk) }] 

. . , -~ (i6) 

.Fig. 4 : is a .diagram, showing the beam forming unit for 
this embodiment of the. present, invention . As can be seen 
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from Fig. 4, the beam forming unit 4 includes an arrival 
direction estimation unit 6, an inverse modulator 8 and a 
beam former 7. The beam former 7 includes multipliers 41 to 
44 'and an adder 45. 
5 The despread output signals x {0) (n) to x l3) (n) 

corresponding to the A antenna elements are input to the 
beam forming unit 4, which are inverse modulated by the 
inverse modulator 8 using- the known transmission symbol such 
as the pilot symbol/ as described above in conjunction with 

10 formula; (10) . As a result, inverse modulated signals zO to 
... z3. are output from the inverse modulator 8. For example, 

the signals zO through z3 can be obtained by multiplying the 
pilot symbol by the despread output signals x <0) (n) to x (3) (n) 
upon reception of the pilot symbol/ so as to obtain signals 

15 corresponding to the carriers of the desired signal 
excluding the pilot symbol component. 

The arrival direction estimation unit 6 calculates 
weighting coefficients aO to a3 based on the inverse 
. .mqduia ted signal^'- \±Q:':'to\ z3V and inputs the weighting 

20 . coefficients aO to; a3 . to .the multipliers 41 to 44 of the 
beam former "7. ; .As'^slipwh in fig. 5, the arrival direction 
estimation unit .6 includes correlation calculators 51 to 53, 
an adcie?: 54 : and an array weighting coefficient calculator 
55 ; ; As shown in Fig-* 6, the correlation calculators 51 to 

25 53 are each provided with a first shift register 61, a 
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second shift register 62, an adder 63, a multiplier 64 and 
"an averaging filter 65. 

The arrival direction estimation unit 6 shown in Fig. 5 . 

has a structure similar to that — — '> shown in 

5. .Fig, However, her e> the correlation calculators are 

constructed differently. In the correlation calculator of 

this embodiment of the present invention, the inverse modulated output signal 

• z- m) (n) corresponding to the mth antenna element is input to 
the first shift register 61, and the inverse modulated 
10 output signal 2 (nrtl) (n) corresponding to the (m+l}th antenna 
■ ./• element is input tb the second shift register 62. The adder 
63 sums symbols L in both the positive and negative 
directions, excluding the symbol corresponding to the zero 
. lag.. from the second. shift register 62. An output of this... 
15 adder 63 is then input to the multiplier 64. 

Accordingly/ the multiplier 64 obtains a complex 
conjugate product of the inverse modulated output signal 
z itn) (n) cor re spending to the mth antenna element obtained via 
the' first -shift .register 61 and the inverse modulated; output 
20 signal z IRH1> (n) corresponding to the (m+1) th antenna element 
for the L symbols in both the positive and. negative 
directions excluding the zero lag at the center of the 
second shift- regist^ output, signal of the 

'^multiplier , 64 is- then- averaged by the. averaging filter .65. 
25 In othier ".Words," it is possible to obtain the correlation R ln>) 
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from the following formula (17). In addition, formula (18) 
is obtained by developing the terms following the summation.: 
Further/ by averaging the correlation R |m) of the formula 
(1?) for each of. the antenna elements, the cross-correlation 
function among the adjacent antenna elements can be 
obtained from. formula (19). 



■." L 

R k «»'«=< Z z lB| (n) z f,?tl) * (n+l)> 
10 ' ■ , . 1-.-L . 

1* 0 (17) 



L 

15 X z ,m ' (n) z ta * 11 * (n+1) 

1=-L 
1*0 . 



20 . .... =lA k | 2 exptj <a+l> )]+Z . . XlAipWi^n) 

' ••• • .. i=-L .- i=0 

1*0 i*k 
25 W llt * (n+l)exp[j ((z>i <mr -0i tB * 1) ) ] 

r'' ■ ,■ : \L '•' : '• • • . ■-. . '• ... 

+ Z N ln) (n)jj 1 ?'* (n+1) — >r (18) 

30 : _ 1*0 :' 

"V- ■-■ . M-2 ' '•' : '" ■ 

> . ; . . . Rk={i/(Wi>..z Rk (n,> — (19) 
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m=0 

in Fig. 5, the output signals R ll0) through of the 
correlation calculators 41 to 53 are then summed by the., 
adder, 54. However, the term [1/ (M-l) ] of formula ( 19) 
averages depending on the number of array antenna elements / 
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which is determined in advance. Accordingly, the function 
of carrying out the process that corresponds to the term 
[1/ (M-l) ] is not included in Fig. 5. 

Using the cross-correlation function R* of the formula 
(19), the arrival angle 6* of the kth user is obtained from 
the following formula (20). 

9 k ^sin^t (l/2*d) tan" a {Im(R K )/Re(R K ) }] 

— (20) 

That is, the arrival angle 9 k is obtained based on the 
ratio of the imaginary part im (R^) and the* real part Re (K k ) 
of the cross^correlation function R*. .7 

The weighting .coefficient a k (n0 input to the beam former 
7 can be obtained from tbe following formula (21) . Further, 
output signals y k , D (n) and y k ,p(n) of the beam former 7 are 
described by the following set of formulas (22) • 

: a Jt (m) -(R x *) ro (m=0 to M-l) : (21) 

y^ofn)^ 2 a k [ ^D k W{n) ";. Data Signal V > 

m=0 

M-l ■ • ' ' ■■ .'':: [■-'. ; V 

:-yk, ? (n)? S,a k <a, x k ,?) .(n) .-: -Pilot -signal-..;'; 

".' ' .m^o : ■ "... 

. ' \ " . : ,. ;,•/:/.} >— i^V 

In other words, in Fig.. 4 the. weighting, coefficients . 
a' 0 ' to' a u > from the arrival direction estimation unit 6 are 
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obtained from formula (21) . The weighting coefficients a {0) 
to a ,3) are respectively input to the multipliers 41 to 44 of 
the beam former 7 and multiplied by the despread output 
signals x l0) (n) through x <3> (n) . The adder 45 then produces an 
output y(n) by summing the products of the inverse despread 
output signals x {0) (n) to x <3} (n) and weighting coefficients 
a I0) to a l3) while matching the phase with the estimated 
arrival direction. The phase is matched with the arrival 
direction because the weighting coefficients a (D) through a <3 
correspond to the complex conjugates of the cross- 
correlation function R*- The output y(n) of the adder 45 is 
then input to the channel receiver 5 (refer to Fig. 1), so 
that : the pilot signal y P and data signal y D corresponding to 
the phannel are demodulated and output . 

It should be noted that the present invention is hot 
limited to the embodiment described above, and various 
variations and modifications may be made. For example, the 
functions of the operation processes of the arrival 
•direction estimation unit 6 and the like can be realized by 
the functions of a digital- signal, processor (DSP) , and other 
functions can similarly be . realized by the functions : of the 
processor . In addition, ., it. is- possible to provide an . 
interference canceller or the like to eliminate the \ 
interference among the channels. 
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?& described abo&, axxxding to an eibodinnent of the present inventioa, 
a D5-CDMA communication system performs; a transmission by a 
spreading process using a code having a period longer than a 
symbol length. Further, received signals from corresponding 
5 array antenna elements 1-0 to 1t-{M-1) are subjected to a 

despreading process." A cross-correlation function is 

also obtained for lags (preferably excluding the zero lag) for signals 
corresponding to adjacent antenna elements obtained by 
inverse-modulating despre ad output: signals by a known 
10 symbol. An arrival direction of a desired signal is 

estimated based on a cross-correlation function. This 
ennbpdiment -his an advantage in that the arrival: direction 

of the desired signal can. be estimated accurately, even in 
; cases where the reception lev^l. of the desired signal is low 
15 or the level of the interference signal is high. Therefore, 
it is possible to increase the antenna gain of the array 
antenna, to eliminate the interference/ and to reduce the 
required transmission power. 

In another eriixdiment , by obtaining the ^cross-correlatiOT function 
. 20 for a; plurality/ o^' s 

directions other : than the zero lag among inverse modulaLted 
signals corresponding to the adjacent antenna elements and 
averaging -the 'crdssrcbrrelation functions, .there -is an 
advantage in that the degree of e's tima t ing * the a:r rival 
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direction of the desired signal can be improved, even with 
respect to a phase deviation caused by fading or the like. 
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CLAIMS : 

. 1. An. arrival direction estimation netted fear use in a E&-am 
communication system . including an array antenna for 
performing a transmission by a spreading process using a . 
5 code . having a period longer than a symbol length, said 
method comprising the steps of: 

subjecting received signals from the array antenna 
.'•• to' a despreading process in order to produce despread output 
: signals; : 

: 10 obtaining a cross-correlation function of time- 

lags of the received signals from adjacent antenna elements, 
.by inverse-modulating the despread output signals by known 
symbols; and 

- - estimating, an arrival direction of one of the. .... 
15. received signals based on the cross-correlation function* 
2. The method of claim 1, wherein said cross* 
. ; correlation function of the time lags between "the received 
. signals, of adjacent antenna elements exclude a zero time 
. ."''■ ! '. ;.'"*,\; ias> \-..V .••'.''.. 

.;;*20 •'/;'. '3 v . The method of claim 1 ctr 2, further comprising the steps 
/'■■•V'V;/ ;;Qf ..,*.'•..',;-■. :.. ,■■ ,-. ."'• v ''-'V'-\ 

• averaging values; of the cross-correlation function 
' for a 'plurality of time lags as the^cfross-correlation 
".. fxinctidh of ■ the time lags of the received signals fi;om the 
25 adjacent .antenna . elements; 
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4. The- arrival direction estimation method of any preceding claim, 
wherein said obtaining the cross-correlation function uses 
both positive and negative time lags between received 
signals 

5. The method of any preceding claim, wherein the known symbol is 
a pilot symbol, 

5„ The method of any one of claims 1 to 4/ further ccnprising the step 

o f : '.' 

providing a decision feedback symbol as the known 
symbol 

7 . An arrival direction estimation method for. use in a DS- 
CDMA communication system using signals sub j ected . to a v 
spreading process by a code, said method comprising the 
steps/of: .■ . 

subjecting the signals to a despreading process, in 
order to produce despread signals; 

inverse modulating the despread signals by a knbwn 
symbol in orcier to obtain inverse modulated signals r and 

producing weighting coefficients corresponding to. 
arrival directions of the signals from the inverse modulated.- 
signals.;' • .\ ' : : ^ ■; 

' . 8. The method of claim 7 , wherein the code has a period 
.Ibhger than a symbol length. 

9. The method of claim 7 or 8, wherein the. producing of 
weighting coefficients includes obtaining a cross- 
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correlation function of time lags from the inverse modulated 
signals and producing the weighting coefficients based ori 
the cross -correlation function. 

10. The method of any one of claims 7 to 9, wherein the cross- 
. 5 correlation of the. time delays excludes a zero time lag. 

11. The method of any one of claims 7 to 10, wherein the knoai syibol is i£ 
a pilot, symbol. . 

12 . The rretird of aiy one of clairrs 7 to 10, vterein the known sytol is 
a decision feedback symbol obtained from a channel receiver. 
10 13. The method of any one of claims 7 to 12, which further includes 

multiplying the despread signals by the weighting 
coefficients to produce weighted signals. 

14. The method of claim 13, which further includes 
. . summing^ the weighted signals in order to produce . output . . 
15 signals. 

15- A DS-CDMA receiver unit for use in a DS-CDMA ccrrrnunication 
• system to produce a transmission by a spreading process, 
using a code having a period longer than a symbol lengthy 
/ said receiver unit comprisingi: 
' '. 20 a de spreading unit for subjecting received signals . 

to a de sp reading process to produde despread 6utput. sigriais; 
■■■■a, beani; f orming. unit for receiving said despread output 
. 1 -signals/ from said despreading unit; anci j 
a ; channel receiver for receiving output signalis v 
25 . from said beam forming unit ; 
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said beam forming unit . including an inverse 
modulator for inverse-modulating the despread output signals 
by a known symbol, an arrival direction estimation unit for 
producing weighting coefficients based on signals. from said 
5. inverse modulator, arid a beam former for producing an output 
by multiplying the weighting coef f icieiits f rom said arrival 
direction estimation unit with the despread output signals 
producing weighted signals and then combining the weighted 
•signals..; 

10 16. The DS-CDMA receiver unit of claim 15, wherein said 

inverse modul ator invers e^modulat es the de spread output 
signals from said despread unit using a, pilot symbol. 

17 1 The ds-CDMA: receiver unit of claim 15, wherein said 
inverse modulator inverse-modulates the despread output 
15 signals from said despreading unit by a decision feedback 
symbol from said channel receiver, 

IB. The DS-CDMA receiver unit of any one of claims 15 to 17, vdnerein said 

arrival direction estimation unit includes a correlation 
calculator for obtaining a cross-correlation of the signals 
20 from said inverse /modulator, and an array weightirigr 
coefficient calculator for calculating the weighting 
coef f icients- based on the cross- co r r ei a t i on obtained by said 
correlation calculator. 
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.19, The DS-CDMA receiver unit of claim 18/ wherein said 
cross-correlation includes a plurality of positive and 
negative time lags excluding a zero time lag, 

2.0. : .The. DSr-CDMA receiver unit of claim 18 or 19, wherein said- 
5. correlation calculator includes a first shift register and a 
second shift register for receiving signals from said 
inverse modulator, a multiplier for obtaining the 
correlation, and an averaging filter, . said second shift 
register provides to said multiplier a plurality of positive 
10 and negative symbols with respect to a symbol input to said 
multiplier from said first shift register. 

21. An arrival direction estimation method substantially as 
hereinbefore described with reference to Figures 1 to 6 of 

the accompanying drawings. . 

22. A CDMA receiver unit substantially as hereinbefore 
described with reference to Figures 1 to 6 of the accompanying 
drawings . 
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